Impurity resonance states in electron-doped high T c superconductors 
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Two scenarios, i.e., the anisotropic s-wave pairing (the s-wave scenario) and the d-wave pairing 
coexisting with antiferromagnetism (the coexisting scenario) have been introduced to understand 
some of seemingly s-wave like behaviors in electron doped cuprates. We considered the electronic 
structure in the presence of a nonmagnetic impurity in the coexistence scenario. We found that even 
if the AF order opens a full gap in quasi-particle excitation spectra, the mid-gap resonant peaks in 
local density of states (LDoS) around an impurity can still be observed in the presence of a d-wave 
pairing gap. The features of the impurity states in the coexisting phase are markedly different from 
the pure AF or pure d-wave pairing phases, showing the unique role of the coexisting AF and d-wave 
pairing orders. On the other hand, it is known that in the pure s-wave case no mid-gap states can 
be induced by a nonmagnetic impurity. Therefore we proposed that the response to a nonmagnetic 
impurity can be used to differentiate the two scenarios. 
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I. INTRODUCTION 

Up to now, the pairing symmetry of electron-doped 
high T c superconductor is still under debate, and vari- 
ous interpretations for experimental results are contro- 
versial. On one hand, the phase sensitive Josephson 
junction experiments indicate d-wave pairing symmetry 
pairing, [1] and meanwhile angle-resolved photoemission 
spectra (ARPES) [2] and electronic Raman spectra[3] 
suggest a non-monotonic d-wave energy gap as a function 
of momentum. On the other hand, tunnelling spectra [4] 
and low temperature specific heat [5] show more or less s- 
wave like behavior, i.e. the absence of zero bias conduc- 
tance peak in the spectrum. Two scenarios have been 
introduced to explain these results. One is the s-wave 
scenario, in which the s-wave character is regarded to be 
intrinsic, but it is incompatible with the phase-sensitive 
Josephson junction experiment. The other one is the 
coexisting scenario in which the coexisting antiferromag- 
netic (AF) order disguises the d-wave character of su- 
perconductivity (SC).[6, 7] The idea is as follows. First 
of all, the in-plane AF order has been observed by neu- 
tron scattering [8] and transport experiments [9], and this 
AF order should be considered naturally. Further, since 
the AF order opens a full gap in single particle excita- 
tion spectrum, which combines with the nodal d-wave 
SC gap to form an effective full single-particle gap. In 
a recent paper it is shown that a more careful analysis 
of the Raman spectra could disentangle the effective gap 
into its AF and pairing components. [10] It remains to see 
whether a single-particle probe, such as the scanning tun- 
nelling microscopy (STM) could be used to identify the 
ingredients of the effective gap and thence differentiate 
the two scenarios. 

In fact, impurity in unconventional superconductors 
proves to be a useful tool to characterize various SC or- 
ders. For example, in a d-wave superconductor a zero 
energy impurity resonant state appears as a hallmark of 



the d-wave pairing symmetry. [11]. On the contrary, in 
conventional s-wave superconductor the resonant states 
lie at the gap edge, which is known as Yu-Shiba-Rusinov 
state. [12] The drastic difference is due to the phase struc- 
ture of the two SC orders: in the d-wave case the phase of 
pairing wave function changes sign across the nodal lines, 
while in s-wave case such no sign change occurs. Mid-gap 
impurity resonant states can also occur in p x + ip y and 
d x + id v superconductors. [13] Impurity resonant states 
can also be used to characterize the electronic structure 
in some materials. [14] Furthermore, it is also proposed 
that by inspecting the line-shape of the resonant peak as 
a function of temperature around a nonmagnetic strong 
impurity can differentiate the phase disorder scenario 
and d-density wave scenario of the pseudo-gap phase in 
hole-doped superconductors. [15] Therefore, STM mea- 
surements of the impurity states can provide important 
messages on the underlying system. [16] The question we 
now ask is if AF and d-wave pairing coexist, what is the 
nature of the impurity state, and in particular whether 
a low energy resonance state can still arise around the 
impurity. 

In this paper we calculate the LDoS around a nonmag- 
netic impurity in electron-doped high T c superconductor. 
Our main results are as follows. Firstly, although the 
AF order gaps the quasi-particles and the bulk density 
of states is s-wave like, two mid-gap impurity resonant 
states, lying symmetrically at positive and negative en- 
ergy, can be observed in the presence of a d-wave SC 
order. Secondly, the two resonant peaks in LDoS ap- 
proach and cross each other when the impurity scattering 
strength increases up to the unitary limit. At an inter- 
mediate scattering strength the two peaks merge into one 
peak at the Fermi energy. Since such a mid-gap impu- 
rity resonant state does not exist in the s-wave scenario, 
the two scenarios for the electron-doped superconductor 
are differentiable by STM measurements of the nature 
of the low energy impurity states. The structure of the 
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rest of the paper is as follows. In section II we describe 
the model and method. In section III we present the re- 
sults. First we obtained the phase diagram of the system 
in the absence of impurities. Second, we discuss analyt- 
ically a toy model with particle-hole symmetry to show 
the existence of impurity resonance states in the coexist- 
ing scenario. Third, we switch back to the actual situa- 
tion in the phase diagram, and calculate the LDoS in the 
presence of a nonmagnetic impurity. In order to see the 
particular role of the coexistence orders, we compare the 
actual case to the cases with AF or d-wave pairing alone. 
Section IV is a summary of the work. 



II. MODEL AND METHOD 

We adopt the t-t'-t"-J model on a square lattice with 
the hamiltonian, 

h = -t (4^ + h - c -) - *' Yl (4^ + h - c -) 



Here Cj CT and c\ a are Fermion operators subject the the 
non-doublc-occupancy constraint, {13)2, an d (ij}3 

denote the first, second and third nearest neighbor pairs 
respectively. It should be noted that we work in the hole 
picture, so that a hole in the above model represents 
a physical electron-double-occupancy in electron doped 
cupratcs. To cope with the above t-J model, we insist 
in calling the electron-double-occupancy as holon unless 
indicated otherwise. For the parameters we choose \t\ as 
the unit of energy, so that t = -1, t' = 0.32, t" = -0.16 
and spin exchange integral J = 0.3. [6] We emphasize that 
the choice of parameters are conventional, but our results 
are not sensitivity to the parameters. 

We apply the slave boson mean field theory (SBMFT), 
within which the projected Fermion operator is decou- 
pled to a spinon and a holon part c iiy = h\fi a , and the 
restriction of no double occupancy is replaced by the con- 
straint h\hi + ^2 a f i(7 jia = 1-[17] In the mean field theory 
this operator constraint is replaced by its average coun- 
terpart. The holons are assumed to condense at zero 
temperature, so that hi,h\ — ► y/x, where x is the doping 
level. The spin-exchange term is decoupled in a stan- 
dard way into hopping, pairing and spin-moment chan- 
nels, [18] 
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-3/8((Aj j )A y - + h.c) + {{rhijrhj + m^rrij)) (2) 

Here the bracket (•) denotes mean value of an opera- 
tor. Xij = J2a fLfj" is tlie no PP m g operator, A^ = 
fafjl — filfjl is the singlet paring operator, and fn t = 



Y] crf} a fi(j is the magnetic moment operator. The mean 
value of these operators are the corresponding order pa- 
rameters. 

We introduce 4-spinors, 

= (/kT,/ikl'/k+QT>/ik-Q|) T 

in the momentum space. Here Q = (it, it) is the nesting 
vector. The mean field Hamiltonian can be written in 
terms of the 4-spinors as H^ MF — 'Jk^k'I'k, where the 
4x4 matrix hk is given by 

I a A k 
A k -a 

2Jm 
2Jm A 



2Jm 

2Jm 

£k+Q A k+ Q 

k+Q -£k+Q j 



(3) 



where £k is the mean-field dispersion of the Fermions, Ak 
is the paring gap function in momentum space, which are 
given by, 

£k = — M - ( 2xt + 3/4x)[cos(fc ;E ) + cos(k y )] 

— Axt' cos (k x ) cos(k y ) — 2xt h '[cos(2fc x ) + cos(2fcj,)], 
A k = -3/4JA[cos(fc 2: ) + cos(fc y )]. 

Here m = (m), \ = (x)i A — (A), and /j, is the chem- 
ical potential to fix the doping level. All of these pa- 
rameters are obtained through self-consistent mean field 
calculation. The spinon's propagator in the above repre- 
sentation is obtained as Gk(iu) n ) = (ico n l — ft-k) -1 - The 
coherent part is given by Q co h — xQ upon convoluting 
with the holon part. 

For later discussion we define a 2 x 2 Green's function 
Goirt^r) = -f((.n h fl) T (r)(f] v f 3l )(0)) for the 2- 

spinors (fa , f^) T in the real space and r is the imaginary 
time. It's Fourier component in frequency uo n is related 
to Q as follows, 
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+Qk,me iQri +g K ive lQ -^- r ^}, (4) 



where the summation is over the magnetic Brillouin zone 
(MBZ) and TV is the size of lattice. The subscripts 
I, II, III, IV indicate the left-top, right-top, left-bottom 
and right-bottom 2x2 block elements of Q. 

The next step is to obtain the Green's function in the 
presence of a single impurity. We model the impurity 
by a single-site potential located at the origin, and adopt 
the T-matrix formulation[13, 15, 20] to calculate the per- 
turbed Green's function G by 

G(r i ,r j ;z) = G (r i ,r j ;z) + G (r i ,0;z)T(z)G (0,r j ;z),(5) 
where 



T(z)= Va 3 + Va 3 G (0,0;z)T{z). 



(6) 
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FIG. 1: (Color online) SBMFT phase diagram of electron 
doped superconductor. The red dash-dotted line, green 
dashed line and blue solid line denote the hopping magni- 
tude, d-wave paring order and staggered magnetic moment 
respectively. The inset shows the bulk density of states at 
x = 0.13. 



Here z is the complex continuation of iu> n , V is the non- 
magnetic impurity potential, and 173 is the third Pauli 
matrix. In the T-matrix formulation the correction of 
mean-field order parameter induced by impurity is ig- 
nored. This was shown to be sufficiently consistent with 
a full self-consistent calculation. [19] 

Finally, the LDoS can be measured by STM directly, 
which is given by 
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FIG. 2: (Color online) The LDoS at position (0,1) for various 
impurity potential strengths in a particle-hole symmetrical 
system. Form bottom to top the potential strengths are suc- 
cessively 0, 5, 10, 15, and 1000, in units of \t\. Vertical offset 
is used for clarity. 



on the coexisting phase. A typical result at doping level 
x = 0.13 is shown as the inset of Fig.l. The U-shaped 
DoS indicates the absence of node in single particle exci- 
tation spectrum, and have been observed in point contact 
tunnelling spectra but was interpreted as the character of 
s-wave pairing. [4] In our case, it is a result of coexisting 
AF and SC orders. Similar behavior was found by Yuan 
et al using t-U-V model. [6] 



B. Impurity states in a particle-hole symmetrical 
case 



where 77 denotes an infinitely small positive number. The 
convention of the sign before to is chosen for the electron- 
doped case under our concern. 



III. RESULTS 
A. SBMFT phase diagram 

We first present the main results of our SBMFT calcu- 
lation. The mean field phase diagram is shown in Fig.l. 
The calculation is done on a 1000 x 1000 lattice. We 
can see that there are two phase transition points. At 
x ~ 0.14 there is a first order phase transition form the 
AF phase to the paramagnetic phase. On the other hand, 
at x 0.20 there is a second-order phase transition from 
the superconducting phase to the normal phase. Thus in 
the range x < 0.14 the AF and SC order coexist. The 
optimal doping level lies at x ~ 0.06. 

Given the mean field order parameter as above, we can 
calculate the bulk density of states (DoS) in different 
phases. In the following discussion we will concentrate 



To see whether resonant impurity states could ap- 
pear in the coexisting phase, we digress to consider a 
the particle-hole symmetrical case that is analytically 
tractable. To reach this case we simply set /.i = 0, t' = 0, 
and t" — 0. We have Ak = — Ak+Q and £k = — £k+Q- 
The unperturbed on-site Green's function Go(0, 0; z) can 
be calculated explicitly as, 



Go(0,0;z) 



1 



(27T) 



dk 1 



BZ 



b 2 



xa (8) 



Here the Co is the two by two unit matrix, and b = 2Jm 
is the Curie- Wess potential. The T-matrix is given by 
T^ 1 {z) = K-^-G^O.O;*). The position of the resonant 
state is given by detT -1 = 0. In the unitary limit V — > 
00 the resonance condition is given by det Go(0, 0; z) = 
0.[13] Since the off-diagonal elements of Go is zero due 
to the d-wave pairing symmetry, the resonance occurs 
at z = b. Since the 11 (22) component of G describes 
particles (holes), we see that another resonance should 
occur at z = —b in the particle picture alone. This is 
also evident from the two contributions in the expression 
for M in Eq.(7). 
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FIG. 3: (Color online) The same plot as Fig.2 but for the 
actual cases described by the phase diagram. 
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FIG. 4: (Color online) The same plot as Fig.2 but in pure AF 
cases at two doping levels. 



If the impurity potential is finite, the situation is more 
complicated. In Fig.2 we present the dependence of LDoS 
on potential strength. In the calculation wc simply set 
X = 0, A = 0.7 and b = 0.3 for illustration. The LDoS 
is for site (0,1), a nearest neighbor of the impurity site. 
When the potential strength is V/t — 5, there appears 
already two resonant peaks lying symmetrically at the 
positive and negative sides already. With increasing po- 
tential strength the resonant peaks shift toward the gap 
edges, and the height of peaks is strongly enhanced. Sim- 
ilar phenomenon has been observed in the study of or- 
ganic superconductivity with bond-dimmerization.[21] In 
the next subsection we discuss how the resonance behaves 
in the realistic situation. 
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C. Impurity states in electron doped cuprates 



FIG. 5: (Color online) The same plot as Fig.2 but in pure d 
wave SC cases at two doping levels. 



In this subsection we discuss the realistic situation in 
electron doped cuprates. We choose the result of SBMFT 
as the input of T-matrix formulation. The evolution of 
LDoS at site (0,1) with impurity potential strength at 
the doping level x — 0.13 is shown in Fig. 3. The mean 
field order parameters are \ — —0.37, A = 0.054, and 
m = 0.19. At first we note that in the unitary limit 
V/\t\ — 1000, there are two mid-gap resonant states lying 
symmetrically at the positive and negative energy, but 
the height of the two peaks are different due to the effect 
of coherence factor in G°(r, 0; z), due to the particle-hole 
asymmetry in the present case. The two peaks are visible 
for V/\t\ > 5. With increasing potential strength, wc 
see that the two peaks seem to cross each other. For 
V/\t\ = 10 the two peaks meet and merge into one peak 
which lies right at the Fermi energy. It is therefore clear 
that resonance states do appear in the coexisting phase, 
even though the nodal d-wave gap is disrupt by the AF 
gap. 



D. Pure AF order and pure SC order case 

In order to understand the unique role of the coexisting 
orders, we compare the results in the pure d-wave case 
and the pure AF order case. We set the d-wave SC order 
parameter to zero in (3) and leave the AF order param- 
eter unchanged for pure AF order case. Conversely, we 
set the AF order to zero for pure d-wave case. In prin- 
ciple one should do the self-consistent calculation again 
in each case. We have done so but found no significant 
changes to the remaining order parameters. 

First, we consider the pure AF case. The result is 
shown in Fig. 4(a) and (b) for x = 0.06 and x = 0.13, 
respectively. For x = 0.06, the AF order is strong enough 
to split the two bands, and impurity resonant peaks lie 
in the gap as well as outside the band. For x — 0.13, the 
AF order is weak and the two AF-order-induced bands 
overlap, and then impurity resonant states (the sharp 
peaks ) lie outside bands. The positions of resonant peaks 
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move toward higher energy for higher potential strengths, 
and are pushed to infinity in the unitary limit. 

Second we turn to the pure d-wave case shown in 
Fig. 5(a) and (b) for x = 0.06 and x = 0.13. We can 
identify two resonant peaks at the two doping levels in 
the unitary limit, but the strong particle-hole asymmetry 
has weaken the strength of peaks. On the other hand, 
the two peaks never cross each other. Moreover, for weak 
potential strength these peaks are indiscernible. 

To close this section, we conclude by comparison that 
two well-defined mid-gap resonant peaks appear and 
cross each other with increasing impurity potential only 
in the coexisting phase with both AF and d-wave SC 
order. 

IV. SUMMARY 

In this work we discuss the impurity resonant state in 
the coexisting phase with d-wave SC and AF orders. We 
demonstrate analytically the existence of mid-gap reso- 
nance states in a putative particle-hole symmetrical case, 
and present the results in the realistic case, where the 
resonance peaks can shift and switch with increasing po- 
tential strength. At an intermediate potential strength, 
the two peaks merge at the Fermi level. These unique 
features do not appear in the separate pure AF or pure 
d-wave pairing states. It is also know that mid-gap impu- 
rity states do not appear in a pure s-wave phase. Thus, 
the impurity state can be used to differentiate the s-wave 
scenario and coexisting scenario. 
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